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A theory is presented which allows us to quantitatively calculate the excess sur- 
face tension of acid solutions. The H + , in the form of hydronium ion, is found to be 
strongly adsorbed to the solution-air interface. To account for the electrostatic poten- 
tial difference measured experimentally, it is necessary to assume that the hydronium 
ion is oriented with its hydrogens pointing into the bulk water. The theory is quan- 
titatively accurate for surface tensions and is qualitative for electrostatic potential 
difference across the air-water interface. 
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I. INTRODUCTION 



Electrolyte solutions are of fundamental interest for a variety of disciplines. Over a hun- 
dred years ago Hofmeister observed a strong dependence of the stability of protein solutions 
on the specific nature of electrolyte. While some ions tend to stabilize protein solutions, 
often denaturing them in the process, others destabilize them favoring protein precipitation. 
A few years after Hofmeister, Heydweiller- observed that salt increases the surface tension 
of the air-water interface. Furthermore, Heydweiller noticed that the relative effect that ions 
have on the surface tension follows closely the Hofmeister series, suggesting that the two 
phenomena are related. 

Over the last hundred years there has been a great effort to understand how ionic speci- 
ficity influences stability of protein solutions and how it affects the surface tension of the 
air-water interface. Langmuir- was probably the first to attempt to construct a quantita- 
tive theory of surface tensions of electrolyte solutions. Appealing to the Gibbs adsorption 
isotherm, Langmuir concluded that the excess surface tension of electrolyte solution was a 
consequence of ionic depletion from the interfacial region. However, no clear explanation 
for this depletion was provided. A few years after, Wagner 3 argued that this depletion was 
the result of interaction between the ions and their electrostatic images across the air-water 
interface. Onsager and Samaras 4 simplified Wagner's theory and obtained a limiting law 
which they argued was universally valid for all electrolytes at sufficiently small concentra- 
tion. More recently, Levin and Flores-MenaA used a direct free energy calculation to obtain 
surface tension of strong electrolytes. To have a quantitative agreement with experiments, 
these authors stressed the fundamental importance of ionic hydration. Nevertheless, the 
theory of Levin and Flores-Mena was not able to predict correctly surface tensions of all 
electrolyte solutions. Bostrom et al. 6 suggested that Hofmeister effect and the ionic speci- 
ficity are a consequence of dispersion forces arising from finite frequency electromagnetic 
fluctuations. This theory predicted that weakly polarizable cations should be adsorbed at 
the air-water interface. This, however, was contradicted by the experimental measurements 
of the electrostatic potential difference^ and by the simulations on small water clusters^"—, 
as well as by the subsequent large scale polarizable force fields simulations^— and the pho- 
toelectron emission experiments^"—. These experiments and simulations showed that some 
anions — and not cations — are present at the solution-air interface^— . To explain this, 
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Levin^ extended the traditional Born theory of ionic solvation to account for ionic polar- 
izability. The new theory predicted that highly polarizable anions can actually prefer an 
interfacial solvation. In a followup work, Levin et al.— >2i used this theory to quantitatively 
calculate the surface tensions and the surface potentials of 10 different electrolyte solutions 
and to reproduce the Lyotropic, Hofmeister, series. 

While almost all salts lead to increase the air-water surface tension, acids tend to lower 
it22i22 The only explanation for this is a strong proton adsorption at the water-air inter- 
face^^— . It is well known that H + ion forms various complexes with water molecules^— . 
It has also been suggested that the high surface adsorption of H + is related to the hydronium 
(H 3 + ) geometry^. This ion, has a trigonal pyramidal structure with the hydrogens located 
at the base of the pyramid™. In this form hydrogens are good hydrogen-bond donors, while 
oxygen is a bad hydrogen-bond receptor—. This favors hydronium ion to be preferentially 
located at the interface, with the hydrogens pointing towards the bulk water and the oxygen 
pointing into the gas phased Explicit solvation energy calculations confirm this picture—. 

The sign of electrostatic surface potential difference is related with the relative popula- 
tion of cations and anions at the interface. Because of high adsorption of hydronium ions, 
one would naturally expect that the electrostatic potential difference across the air-water 
interface for acid solutions should be positive. The experiments, however, show that the 
surface potential difference for acids has the same sign as for halide salts, i.e. is predom- 
inantly negative 7,23 . Frumkin 7 suggested that this apparently strange behavior might be 
a consequence of the incomplete dissociation of acid molecules. A different explanation 
was advanced by Randies— who argued that presence of hydroniums at the interface leads 
to a preferential orientation of water molecules resulting in a dipole layer with a negative 
electrostatic potential difference across it. This conclusion is in agreement with the theory 
proposed in the present paper, as well as with the recent molecular dynamics simulations—. 

In this paper we present a theory that allows us to quantitatively calculate surface ten- 
sions of acid solutions using only one adjustable parameter related to the strength of the 
hydronium adsorption to the interface. Predictions of the theory are compared with the 
experimental measurements. The theory is then used to estimate the electrostatic potential 
difference across the water-air interface for various acid solutions. 
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II. MODEL AND THEORY 



We consider an acid solution in a form of a drop of radius R, where r = R is the position 
of the Gibbs dividing surface (GDS) 20 ' 34 . The water and air will be modeled as uniform 
dielectrics of permittivities e w = 80 and e D = 1, respectively. The surface tension can be 
obtained by integrating the Gibbs adsorption isotherm equation: 

d7 = —T + dp + — T_d/i_ , (1) 

where p± = fceTln^A^) are the chemical potentials and A± are the de Broglie thermal 
wavelengths. In this equation + sign corresponds to the hydronium ion, and — sign to the 
anion. The bulk ion concentration is q, = p+(0) = P-(0), where p±(r) are the ionic density 
profiles. The ion excess per unit area due to existence of the interface is 
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where N is the total number of acid "molecules". The ionic density profiles, p±(r), will be 
calculated using a modified Poisson-Boltzmann (mPB) equation, as discussed later in the 
paper. 

Anions are divided into two categories: kosmotropes and chaotropes. The theory of elec- 
trolyte solutions^ showed that chaotropes, Br~, I~, NO J, and C107 , lose their hydration 
sheath near the GDS and are partially adsorbed to the interface. On the other hand kos- 
motropes, F~, CI - , and SO 2- remain hydrated in the interfacial region and are repelled 
from the GDS. 

To bring an ion of radius to distance z > ah from the GDS requires 5 : 

W (z; a h ) = I" dkc~ 2s(z - a » ) k[sCOSh ( kah) ~ ksinh ( ka h)] ) (3) 
2e w J s[s cosh(kah) + k sinh(A;o/ l )] 

of work. In this equation s = \J (k 2 + k 2 ) and k = ^J^Tiq 2 ci ) /e w kBT is the inverse De- 
bye length. The kosmotropic ions remain strongly hydrated in the interfacial region and 
encounter a hardcore-like repulsions from the GDS at a distance of one hydrated ionic ra- 
dius. On the other hand, strongly polarizable chaotropic anions (Br - , I - , NOj and ClOj) 
loose their hydration sheath and can move cross the water-air interface. However, to avoid 
the large electrostatic energy penalty of exposing the charge to a low-dielectric (air) en- 
vironment, the electronic charge density of a chaotropic anion redistributes itself so as to 



remain largely hydratecU^. The fraction of ionic charge which remains inside the aqueous 
environment, x(z), can be calculated by minimization the polarization energy^ 
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9{z) [7r-9(z)]e J 
In the above equation a is the relative polarizability defined as a = 7i/ajj, where 7« is the ionic 
polarizability, a is the unhydrated (bare) radius, and 9(z) = arccos[— z/ao\. Performing the 
minimization we obtain 
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where Xb = <f /twksT is the Bjerrum length. 

The force that drives chaotropic ions towards the interface results from water cavitation. 
To introduce an ion into an aqueous environment requires creating a cavity, which perturbes 
the hydrogen bond network of water molecules. For small ions, the free energy cost of forming 
a cavity is entropic and is proportional to the volume of the void formecP^. As the ion 
moves across the GDS, its cavitational free energy decreases. This results in a short-range 
attractive potential between the anion and the GDS 20 : 

ua^ for z > a , 

2 / \ (6) 
[2--) for - a < z < a , 

where v ~ 0.3/csT/ A 3 is obtained from bulk simulations^.. For hard (weakly polarizable 
ions) the cavitational free energy gain is completely overwhelmed by the electrostatic free en- 
ergy penalty of moving ionic charge into the low dielectric environment. For soft polarizable 
ions, however, the electrostatic penalty is small, since most of the ionic charge remain inside 
the aqueous environment. The total potential of a soft anion, therefore, has a minimum in 
the vicinity of the GDS, see Fig. [TJ 

The H + ions (protons) do not exist as a separate specie in water. Instead they form 
complexes with water molecules, HsO" 1- and H20^"— ~— . Because of its favorable geometry 
(trigonal pyramidal), the hydronium ion (H30 + ) adsorbs to the water-air interface with a 
preferential orientation^ of oxygen towards the air. We model this attraction by a square 
well potential with a range of a hydrogen bond 1.97 A, 



hyd { 



for z > 1.97 A , 

-3.05 k B T for < z < 1.97 A . 



(7) 




FIG. 1. Potentials for all ions at 1M. The GDS is at z = A. The bulk cavitational potential is 
not considered for kosmotropes, since it does not change along the drop. 

The depth of the potential is then adjusted to obtain the experimentally measured surface 
tension of HC1. The same potential is then used to calculate the surface tensions of all 
other acids. We should stress, however, that one should not attach too much meaning to 
the specific value of the potential depth. The real proton transfer is a quantum mechanical 
process, therefore there is bound to be some arbitrariness in how one models it at a classical 
level. Here we have chosen the range of the square well potential to be one hydrogen bond. 
If one changes this distance, the depth of the potential will have to be modified to obtain 
an optimal fit of the surface tension of HC1 solution. However once this is done, the values 
of the surface tension of other acids will not be significantly affected. Thus, the strength of 
H + potential is the only free parameter of the theory. The total potential felt by H + is then, 
U H {z) = U hyd (z) + W(z; 0), see Fig. HJ 

While the kosmotropic anions feel only the potential W(z; a^) and the hardcore repulsion 
from the GDS, the chaotropic anions are influenced by the total potential^: 



U tot (z) 



W(z; oq) + val + 



for z > do 



W(a ; a )z/a + U p (z) + U cav (z) for < z < a 
U p (z) + U cav (z) for - a < z < . 



(8) 



In Fig. [TJ we plot the potentials felt by various ions at 1M concentration, as a function of 
the distance from the GDS. 
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TABLE I. Ion classification into chaotropes (c) and kosmotropes (k). Effective radii (hydrated or 
partially hydrated) for kosmotropes and (bare) for chaotropes, for which we have also include the 
polarizabilities from Ref. 38 . 



Ions 


chao/kosmo 


, 9 . 

radius (A) 


polarizability (A 3 ) 




k 


3.54 




cr 


k 


2 
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Br~ 


c 


2.05 


5.07 


r 


c 


2.26 


7.4 


no 3 - 


c 


1.98 


4.48 


cioj 


c 


2.83 


5.45 


so 2 A - 


k 


3.79 


* 



The ionic density profiles can now be obtained by integrating the mPB equation: 
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where O is the Heaviside step function, pl hao (r) is the density profile for chaotropic anions 
and p^ os (r) for kosmotropic ones. 

Once the ionic density profiles are calculated, the surface tensions can be obtained by 
integrating the Gibbs adsorption isotherm (eqUJ. The ionic radii and polarizabilities are 
the same as were used in our previous work on surface tension of electrolyte solutions^^i. 
In Tabled we summarize this data. 

The depth of the potential Uh y d(z) (eqEJ) is adjusted to fit the HC1 experimental data^ 
(see Fig. this is the only adjustable parameter of the theory. We find that a square 
well potential of depth —3.05 k^T results in an excellent fit of the experimental data for 
HC1 in the range of concentrations from to 1 M. The excess surface tension of all other 
acids is then calculated using the same potential Uh y d{z). The predictions for the surface 
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FIG. 2. Surface tensions for HF, HC1, HBr, and HI. The symbols are the experimental data for 
HC1— and the lines are the results of the present theory. 

tensions of HF, HBr, and HI are plotted in Fig. |2j Unfortunately, we have no experimental 
data to compare for these halogen acids. For H2SO4 and HNO3 (see Fig. [3]), we find a 
good agreement between the theory and experiment. For HCIO4 the theory overestimates 
the surface tension. This is similar to what was found for sodium perchlorate salt 21 . The 
difficulty is that ClOj is a large weakly hydrated ion. Since the cavitational energy grows 
with the cube of ionic radius, a small error in radius leads to a big error in surface tension. 
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FIG. 3. Surface tensions for H2SO4, HNO3, and HCIO4. The symbols are the experimental data 22 
and the lines are the results of the present theory. 



Finally, we use the theory to calculate the electrostatic potential difference across the 
solution-air interface. The surface potential difference A% = (p(R + a ) — 0(0) predicted 
by the present theory has a wrong sign compared to Frumkin's experimental measurements 
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TABLE II. Surface potential differences for various acids at 1M concentration. Contributions from 
electrolyte and aligned water dipoles. 



Acids 


calculated [mV] 


Frumkin 7 [mV] 


HF 


85.5 


-71 


HC1 


1.24 


-23 


HBr 


-95 


-34 


HI 


-144.8 


-61 


HN0 3 


-84.4 


-48 


HCIO4 


-412 


-82 



- positive instead of negative 7 " 1 ' 23 . Positive sign reflects a strong adsorption of hydronium 
ions to the GDS. The simple dielectric continuum theory presented here, however, does not 
account for the structure of the interfacial water layer. Since the hydronium ion at the 
GDS has a preferential orientation with the hydrogens pointing towards the bulk, presence 
of many such ions will result in a dipole layer. Note that in the absence of hydroniums, 
the water dipoles predominantly point along the interface 39 . The hydronium layer produces 
an electric field E = AirpNhj 'e dA, where Nh is the number of hydroniums at the interface, 
p is the water dipole moment, d is the dipole length, and A is the interfacial area. If we 
suppose that all the hydroniums are perfectly aligned, the potential difference across the 
dipolar layer will be Ax w — — 47rpr + /e . Using the dipole moment of a water molecule, 
p = 1.85 D, we obtain the dipole layer contribution to the overall electrostatic potential 
difference. Adding this to A%, we obtain the total electrostatic surface potential difference 
across the solution-air interface. In Table [III we list the surface potentials of various acids 
at 1M concentration. Clearly these values are an exaggeration of the total electrostatic 
potential difference across the interface, since at finite temperature there will not be perfect 
alignment of interfacial hydronium ions. Nevertheless, the theory should provide us an order 
of magnitude estimate of the electrostatic potential difference. In fact, for most acids we find 
a reasonable agreement between the predictions of the theory and Frumkin's experimental 
measurements 7 . A noticeable exception is the HF. Experimental potential for hydrogen 
fluoride measured by Frumkin is negative, while we find a large positive value. Frumkin's 
value for HF is clearly outside the general trend for halogen acids. In his classical review 
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of electrolyte solutions, Randies^ did not mention Frumkin's result for HF acid, while 
discussing his other measurements. We can only suppose that Randies also did not have a 
complete confidence in this particular value. Experimental measurements of excess surface 
potentials are very difficult. This is probably the reason why Frumkin's measurements of 
surface potentials of acids have not been repeated in over 90 years. 

III. CONCLUSIONS 

In this paper we have developed a theory for surface tensions of acid solutions. The 
hydronium adsorption to the interface was modeled by a square well potential, the depth 
of which is the only adjustable parameter of the theory. The agreement between the theory 
and experiments is very reasonable for different acid solutions at concentrations from 
to 1M. In order to account for the experimental values of the excess electrostatic surface 
potential, we must require a preferential orientation of hydronium ion at the interface, with 
the hydrogens pointing into the bulk. With this assumption we get a qualitative agreement 
with the experimental measurements of the excess electrostatic potentials of various acid 
solutions. At the moment this is the only theory that can account (quantitatively) for the 
surface tensions and (qualitatively) for the surface potentials of acid solutions. 

IV. ACKNOWLEDGMENTS 

This work was partially supported by the CNPq, INCT-FCx, and by the US-AFOSR 
under the grant FA9550-09- 1-0283. 

REFERENCES 

1 A. Heydweiller, Ann. Phys. (Leipzig) 33, 145 (1910). 
2 I. Langmuir, J. Am. Chem. Soc. 39, 1848 (1917). 
3 C. Wagner, Phys. Z. 25, 474 (1924). 

4 L. Onsager and N. N. T. Samaras, J. Chem. Phys. 2, 528 (1934). 

5 Y. Levin and J. E. Flores-Mena, Europhys. Lett. 56, 187 (2001). 

6 M. Bostrom, D. R. M. Williams, and B. W. Ninham, Langmuir 17, 4475 (2001). 

7 A. Frumkin, Z. Phys. Chem. 109, 34 (1924). 

10 



8 N. L. Jarvis and M. A. Scheiman, J. Phys. Chem. 72, 74 (1968). 
9 L. Perera and M. L. Berkowitz, J. Chem. Phys. 95, 1954 (1991). 
10 L. X. Dang and D. E. Smith, J. Chem. Phys. 99, 6950 (1993). 
n S. J. Stuart and B. J. Berne, J. Phys. Chem. A 103, 10300 (1999). 
12 P. Jungwirth and D. J. Tobias, Chem. Rev. 106, 1259 (2006). 
13 P. Jungwirth and D. J. Tobias, J. Phys. Chem. B 106, 6361 (2002). 
14 D. Horinek, A. Herz, L. Vrbka, F. Sedlmeier, S. I. Mamatkulov, and R. R. Netz, Chem. 

Phys. Lett. 479, 173 (2009). 
15 M. A. Brown, R. D'Auria, I. F. W. Kuo, M. J. Krisch, D. E. Starr, H. Bluhm, D. J. Tobias, 

and J. C. Hemminger, Phys. Chem. Chem. Phys. 10, 4778 (2008). 
16 G. Markovich, S. Pollack, R. Giniger, and O. Cheshnovski, J. Chem. Phys. 95, 9416 (1991). 
17 S. Ghosal, J. Hemminger, H. Bluhm, B. Mun, E. L. D. Hebenstreit, G. Ketteler, D. F. 

Ogletree, F. G. Requejo, and M. Salmeron, Science 307, 563 (2005). 
18 B. Garrett, Science 303, 1146 (2004). 
19 Y. Levin, Phys. Rev. Lett. 102, 147803 (2009). 

20 Y. Levin, A. P. dos Santos, and A. Diehl, Phys. Rev. Lett. 103, 257802 (2009). 
21 A. P. dos Santos, A. Diehl, and Y. Levin, Langmuir 26, 10785 (2010). 
22 P. K. Weissenborn and R. J. Pugh, J. Colloid Interface Sci. 184, 550 (1996). 
23 J. E. B. Randies and D. J. Schiffrin, Trans. Faraday Soc. 62, 2403 (1966). 
24 L. M. Levering, M. R. Sierra-Hernandez, and H. C. Allen, J. Phys. Chem. C 111, 8814 
(2007). 

25 M. Mucha, T. Frigato, L. M. Levering, H. C. Allen, D. J. Tobias, L. X. Dang, and P. Jung- 
wirth, J. Phys. Chem. B 109, 7617 (2005). 

26 M. K. Petersen, S. S. Iyengar, T. J. F. Day, and G. A. Voth, J. Phys. Chem. B 108, 14804 
(2004). 

27 P. Petersen and R. J. Saykally, J. Phys. Chem. B. 109, 7976 (2005). 

28 S. S. Iyengar, T. J. F. Day, and G. A. Voth, Int. J. Mass Spectrometry 241, 197 (2005). 

29 M. Eigen, Angew. Chem., Int. Ed. Engl. 3, 1 (1964). 

30 G. Zundel, Adv. Chem. Phys. Ill, 1 (2000). 

31 D. Marx, M. E. Tuckerman, J. Hutter, and M. Parrinello, Nature 397, 601 (1999). 
32 L. X. Dang, J. Chem. Phys. 119, 6351 (2003). 



11 



Randies, J. E. B., Advances in Electrochemistry and Electrochemical Engineering, vol. 3 
(Interscience, New York, 1963). 

C. H. Ho, H. K. Tsao, and Y. J. Sheng, J. Chem. Phys. 119, 2369 (2003). 
K. Lum, D. Chandler, and J. D. Weeks, J. Phys. Chem. B 103, 4570 (1999). 

D. Chandler, Nature (London) 437, 640 (2005). 

S. Rajamani, T. M. Truskett, and S. Garde, Proc. Natl. Acad. Sci. U.S.A. 102, 9475 
(2005). 

N. C. Pyper, C. G. Pike, and P. P. Edwards, Mol. Phys. 76, 353 (1992). 

S. M. Kathmann, I. F. W. Kuo, and C. J. Mundy, J. Am. Chem. Soc. 130, 16556 (2008). 



12 



